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We  have  measured  the  x  and  >•  components  of  both  the  intrinsic  and  induced  second-order  nonlinear  susceptibilities  in  glass 
optical  fibers,  in  addition.  v,e  have  measured  the  angular  dependence  of  the  intrinsic  and  induced  second-harmonic  light  intensity 
as  a  function  of  input  polarization.  Our  results  are  discussed  in  relation  to  the  models  put  forward  to  date  to  explain  second- 
harmonic  generation  in  optical  fibers. 


1.  Introduction 


To  date,  several  models  have  been  put  forward  [  1- 
5], -in  an  attempt  to  explain  second-harmonic  gen¬ 
eration  (SHG)  in  glass  optical  fibers..  Due  to  the 
complexity  of  the  glass  structure  most  of  these  models 
are  based  upon  fairly  simplistic  assumptions.  We 
have  in  this  paper  measured  some  fundamental 
properties  of  the  photoinduced  SH  light,  in  order  to 
further  elucidate  the  mechanism  behind  SHG  in  glass 
fibers. 

The  fundamental  properties  we  have  chosen  to 
study  are  the  relative  magnitudes  of  the  tensor  com¬ 
ponents  of  the  intrinsic  and  induced  second-order 
nonlinear  susceptibilities.  The  intrinsic  SH  light  is 
from  electric  quadrupole  interactions  at  the  core¬ 
cladding  interface  of  the  fiber  as  predicted  by  Ter- 
hune  and  Weinberger  [6]  and  recently  verified  ex¬ 
perimentally  [7],  The  electric  quadrupole  tensor, 
which  is  of  rank  four,  has  only  three  independent 


1  Permanent  address:  Thayer  School  of  Engineering,  Dart¬ 
mouth  College.  Hanover,  NH  03755.  USA. 


components  in  glass;  due  to  the  constraints  imposed 
by  the  waveguiding  properties  of  the  fiber  only  the 
ratio  between  two  of  them  can  be  measured  in  this 
geometry.  Measurements  of  the  induced  susceptibil¬ 
ity  are  also  limited  by  the  geometry  of  the  wave¬ 
guide,  since  the  overlap  integral  is  zero  for  certain 
polarization  combinations. 

In  addition  we  have  measured  the  angular  depen¬ 
dence  of  the  intrinsic  and  induced  SH  light  intensity 
as  a  function  of  input  polarization.  The  data  suggest 
that  the  induced  SH  light  comes  primarily  from  the 
initial  electric  quadrupole  interaction  and  that  growth 
occurs  via  improved  phasematching.  It  has  already 
been  suggested  [8]  that  an  increase  in  coherence 
length  could  explain  the  increase  in  the  intensity  of 
the  SH  light.  The  justification  of  such  a  statement  is 
easily  obtained  from  an  order  of  magnitude  calcu¬ 
lation.  In  the  literaiure  it  has  been  reported  that  the 
coherence  length  can  be  increased  from  30  pm  to  30 
cm  [9]  in  specially  prepared  fibers.  This  in  itself 
would  guarantee  an  increase  of  the  SH  light  by  S  or¬ 
ders  of  magnitude  -  in  excellent  agreement  with 
experiments. 


0030-4018/9I/S03.50  ®  1991  -  Elsevier  Science  Publishers  8.  V.  (  Nonli-Holbnd  ) 


549 


9  1 


his  d  acuities.:  ho 3  boon  'approved 
lor  public  reiee.so  end  sole,  ita 
distribution  is  uniunitsdl 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


Volume  S2,  number  5,6 


OPTICS  COMMUNICATIONS 


I  May  1991 


2  Theory 

One  of  the  problems  in  measuring  a  nonlinear  sus¬ 
ceptibility  in  a  waveguide  is  the  interaction  between 
different  modes.  Experimentally,  one  can  measure 
an  effective  nonlinear  susceptibility  which  con¬ 
sists  of  the  susceptibility  times  a  factor.  /  known  as 
the  overlap  integral.  A  good  review  of  overlap  in¬ 
tegrals  in  waveguides  can  be  I'ourd  in  ref.  [  10}.  Fer 
the  modes  in  an  optical  fiber  the  effective  ^jy’  can 
be  expressed  as 

~  J  F»X  Hz,*  d.v  dr  ,  (1) 

Ax. 

where  dardv  is  a  plane  perpendicular  to  the  direction 
of  propagation  and  .V  is  the  normalization  factor.  The 
polarization  P  can  be  written  as 

P.  -  E:Ex  C-) 

for  an  electric-dipole  interaction  and  as 

P,=r$E,VkE,  (3) 

ldr  an  electric  quadrupole  interaction.  In  our  cal¬ 
culations  we  assume  the  nonlinear  susceptibility  is 
constant  and  therefore  move  it  outside  the  integral. 
The  remaining  integral  is  the  overlap  integral  de¬ 
noted  f  which  we  separate  into  a  bulk  part  and  an 
interface  part.  The  bulk  part  corresponds  to  the  core 
of  the  fiber  and  the  interface  part  is  the  region  be¬ 
tween  the  core  and  the  cladding  estimated  from  the 
Goos-Hanschen  shift.  It  is  important  to  note  that  the 
interface  and  bulk  susceptibilities  have  different  units 
related  through  the  thickness  of  the  interface  region 
[  1  1 }.  In  tables  1  and  2  we  have  listed  the  calculated 


Table  I 

Calculated  overlap  integrals  for  electric  dipole  inter?  "lions 


■X-polarization 

polarization 

(v-  HE, 

ut- He,,:  2ru-HEn  core 

E,F. , 

3.03 

3.03 

E.E. 

-6.06X10-’ 

6.06  X  10'3 

(it-  H  E*  1  1 

.  oa-He, 2</>He;,  inlerface 

EJ T, 

0.54 

0.54 

ExEy 

-S.90X  lO'4 

-8.90x  I0'4 

E.E. 

-  I.05x  I0_! 

-3.5Sx  10-’ 

Table  2 

Calculated  overlap 

integrals  for  electric 

quadrupole  interactions 

.x-poiarization 

V’-polarization 

a>Hei 

-VHe;i  core 

F.,E, 

-4.67 

_  i 

EyF.x 

-1.50X  lO'3 

-1.50X  10-’ 

E.E, 

9.18X  10s 

0 

Interface 

E,E, 

-0.24 

-8.07x10--' 

E.E, 

5.ISX  I0-4 

5. !  S  x  1 Q-4 

E.E , 

2.06  X  10! 

0 

values  of  the  overlap  integrals  for  electric  dipole  and 
electric  quadrupole  interactions  separately.  We  have 
limited  ourselves  to  LP0,  (HE,,)  and  LP,,  (HE;i, 
TMji  and  TE0) )  modes  since  they  were  the  only  ones 
observed  in  our  fibers.  The  overlap  integrals  for  po¬ 
larization  combinations  not  listed  in  tables  1  and  2 
are  zero.  In  table  2,  although  'he  electric  quadrupole 
interactions  are  described  via  a  fourth  rank  tensor 
(in  principle  requiring  four  subindexes)  we  have 
limited  ourselves  to  three  subindexes  in  the  table  for 
practical  reasons:  due  to  the  symmetry  constraints  of 
the  tensor,  such  a  reduction  in  indexes  is  unambig¬ 
uous.  Some  of  these  overlap  integrals  (the  ones  re¬ 
lated  to  £(0/0c)£  terms)  might  seem  to  have  un¬ 
usually  large  magnitudes.  In  these  cases,  the  :- 
derivative  multiplies  the  “geometrical"  overlap  in¬ 
tegral  with  the  propagation  constant 

Another  problem  of  measuring  tensor  susceptibil¬ 
ities  in  optical  fibers  is  that  in  a  circular  fiber  every 
mode  of  the  light  has  three  directions  of  polarization 
(£,.  £,.  and  £.).  This  means  that  even  though  our 
input  light  is  linearly  polarized,  in  the  output  we 
inevitably  measure  a  linear  combination  of  tensor 
components. 

There  are  two  ways  to  simplify  the  situation:  one 
is  to  discard  all  the  terms  w  ith  small  values  of  the 
overlap  integral,  and  the  second  is  to  use  a  polari¬ 
zation  preserving  fiber.  From  the  calculations  of 
overlap  integrals,  it  was  obvious  what  terms  could  be 
discarded:  to  implement  the  second  option  we  used 
fibers  with  extinction  ratios  between  the  Ex  and  F, 
polarizations  that  varied  between  0. 1  and  <  10~5.  In 
this  wav  we  could  simultaneously  study  the  influence 
of  the  polarization  properties  of  the  preparation  beam 
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on  the  SH  light  coming  out  of  the  fiber.  However,  it 
is  important  to  remember  that  even  in  the  case  of  a 
polarization  preserving  fiber  (Ey/Ex  <  10'3).  the  ra¬ 
tio  E./F.x  is  approximately  the  same:  for  our  fibers, 
E. <  I0~5£t.  This  meant  that  we  could  simplify  the 
final  expression  to  read  7T<y  =  A:-(r )*^iL 7 £ 
and  similarly  for  l-_u>. 

If  one  attempts  to  measure  an  off-diagonal  tensor 
element  by  exciting  the  fiber  w  ith  light  polarized  45 ; 
to  the  A'-direotion  it  becomes  more  cumbersome.  If 
we  go  through  the  equivalent  derivations  for  this  case, 
it  turns  out  that  the  only  contributing  factor  to  /'(u 
is  the  same  as  for  the  equation  above  but  with 
In  other  words,  because  of  zero-,  .lued 
overlap  integrals,  we  cannot  measure  tensor  ele¬ 
ments  such  as  I  'xWy  in  an  optical  fiber.  To  measure 
r'y.’..  (and  similir  off-diagonal  tensor  elements)  we 
would  have  to  perform  our  measurements  in  a  bulk 
glass  sample.  For  electronic  dipole  interactions  we 
have  similar  difficulties  with  overlap  integrals  being 
zero.  This  is  why  we  have  omitted  the  x.xx/xxy...  etc. 
ratios  in  table  4. 

It  should  be  pointed  out  that  from  these  measure¬ 
ments  alone,  it  is  in  principle  impossible  to  separate 
one  of  the  components  of  the  forbidden  core  second 
harmonic  generation  from  the  interface  contribution 
[  12],  Our  rationale  for  discussing  the  core  and  in¬ 
terlace  contributions  separately  are:  (i),  we  have 
previously  given  strong  evidence  that  the  intrinsic 
SH  light  is  from  a  nonlinear  interaction  at  the  core¬ 
cladding  interface  ( 7  ] ,  and  ( ii )  we  have  also  shown 
that  the  magnitude  of  the  intrinsic  SH  light  was  pro¬ 
portional  to  the  magnitude  of  the  relative  refractive 
index  difference  at  the  interface  favoring  an  electric 
qundrupole  interaction.  Also,  from  a  strictly  math¬ 
ematical  point  of  view  [  12],  it  is  not  clear  that  it  is 
meaningful  to  separate  the  interface  contributions 
from  dipole  and  quadrupole  interactions  as  we  have 
done  in  tables  1  and  2.  However,  following  the  anal¬ 
ysis  ot  Guyot-Sionnest  et  al.  [13],  we  believe  that 
separating  the  interface  overlap  integrals  into  dipole 
and  quadrupole  interactions  gives  more  physical  in¬ 
sight  to  what  is  going  on  at  the  interface. 


3.  Experiment 

Our  measurements  were  done  with  a  (9-switched 
and  modelocked  Nd :  YAG  laser  at  1.2  kHz  and  76 


MHz  respectively.  The  length  of  the  fibers  measured 
was  approximately  7  cm. -We  used  two  different  types 
of  fiber,  one  with  V,„=  2.49  and  V2, „  =  5. 1 6  and  a  sec¬ 
ond  fiber  with  '  3.50  and  7.20,  where  I' is 

the  normalized  frequency.  The  input  polarization 
could  be  arbitrarily  varied  using  a  polarizer  and  a  /./ 
2  waveplate.  For  each  fiber  we  characterized  the  lin¬ 
ear  polarization  properties  both  before  and  after 
“preparation”.  The  polarization  extinction  ratio  var¬ 
ied  between  0.1  and  <  10-3  for  the  different  pieces 
of  fiber  used  at  both  the  fundamental  and  SH  wave¬ 
lengths.  The  preparation  was  done  at  many  different 
laser  intensities  ranging  from  0.5  to  20  k\V.  All  fibers 
w-ere  read  out  at  very  low  fundamental  input  powers 
to  ensure  that  no  growth  took  place  during  the  mea¬ 
surements  Although  the  fibers  were  always  satu¬ 
rated  in  the  A'-direction  (preparation  polarization  for 
all  fibers)  the  read-out  was  made  with  the  input  po¬ 
larization  in  a  new  direction  ( in  which  the  fiber  could 
be  prepared).  The  various  tensor  components  were 
measured  by  setting  the  input  (read-out)  polariza¬ 
tion  to  0  =  ,  45"  and  90:,  while  for  each  of  these  input 
polarizations  the  output  signal  was  recorded  with  the 
analyzer  at  0;  or  90 :.  (The  .v-direction.  preparation- 
direction.  corresponds  to  0:.)  The  result  of  the  dif¬ 
ferent  intensity  ratios  from  these  measurements  can 
be  found  in  table  3  for  the  intrinsic  and  photoin- 
duced  SH  light.  In  table  4  we  have  listed  the  ratios 
between  the  equivalent  tensor  components  that  we 
could  calculate  from  tables  1,  2  and  3.  The  calcu¬ 
lations  were  done  using  the  following  formula: 


x'dl  f*  Ki'/rz)  \fr) 


(4) 


where  /’corresponds  to  an  overlap  integral  and  /". 


Table  3 

Measured  intensity  ratios  for  intrinsic  and  indirect  SH  light  after 
the  optical  fiber 


Ratio  of 

measured 

intensities 

Intrinsic 

Induced 

2.4;  1.2 

30;  10 

2.4±  1.2 

1 2  +  3 

2.7;  1 

3;  1 

2.7  y  i 

8  ;  2 

M,/G, 

1.5  i  0.3 

1 55  i  57 
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Table  4 

Measured  ranos  between  various  components  ol‘ intrinsic  and  in¬ 
duced  second-order  nonlinear  susceptibilities 


Ratio  of 

tensor 

elements 

Intrinsic 

Induced 

- 

XXX  /  I’XV 

1.5  aC. 4 

5.5;  1. 0 

xxx/xvv 

1.5r0.4 

3.5;0.5 

xxx/yyy 

1.2  ±0.3 

1 2  t2 

/:'u  arc  the  measured  light  intensities  for  the  fun¬ 
damental  and  second-harmonic  light,  respectively. 

We  also  measured  the  angular  dependence  of  the 
intrinsic  and  induced  SH  light  as  a  function  of  input 
1R  polarization  as  well  as  output  SH  polarization. 
These  measurements  were  performed  in  the  follow¬ 
ing  way.  To  measure  the  angular  dependency  of  the 
SH  light  as  a  function  of  the  input  polarization  we 
kept  the  SH  analyzer  after  the  fiber  in  the  .v-direction 
(6=  0°)  and  varied  the  input  IR  polarization  in  in¬ 
crements  of  10%  In  this  measurement  r2w  = 
/;cu(03 )  /  /;tu( 90' )  10  for  the  induced  light.  The 

second  kind  of  measurement  was  to  keep  the  input 
polarization  fixed  in  the  .v-direction  and  vary  the  SH 
analyzer  after  the  fiber  in  increments  of  10%  with 
the  result  that  r2oj ~  30  for  the  induced  light.  The  an¬ 
gular  dependence  for  the  two  measurements  is  sum¬ 
marized  in  eq.  (5), 

1(f))  =s/min[  1  +  r:<1,cosJ(d>)]  .  (5) 

The  variation  of  the  intrinsic  SH  light  for  different 
input  polarizations  also  has  a  cos 4(0)  dependence 
on  the  input  polarization  with  r?:2.4.  We  would  like 
to  stress  the  fact  that  the  ratios  given  in  table  3  were 
independent  of  the  preparation  intensity  (the  inten¬ 
sity  was  varied  by  an  order  of  magnitude)  and  the 
degree  of  polarization  preservation  ( the  polarization 
preservation  varied  by  two  orders  of  magnitude)  of 
the  fiber.  Another  importin'  parameter  in  our  mea¬ 
surements  is  the  number  of  modes  at  the  fundamen¬ 
tal  wavelength  that  propagate  through  the  fiber,  since 
our  calculations  are  contingent  upon  the  fundamen¬ 
tal  light  being  in  the  HE,,  mode.  From  measure¬ 
ments  of  the  far-field  angle  we  estimate  that  less  than 
I0°'o  of  the  fundamental  light  could  be  in  a  higher 
order  mode.  This  check  should  always  be  made  when 
short  pieces  of  fiber  (<  1  m)  are  being  used.  The 


reason  is  that  a  fiber  which  is  nominally  ‘‘single- 
mode"  will  slill-transmit  ~80°'o  of  the  light  cxci’ed 
in  the  LP,,  mode  if  the  fiber  is,  e.g.,  ~  10  cm  long 
[14]. 

4.  Discussion 

AH  models  fui  »<iiu  to  date  to  explain  second- 
harmonic  generation  in  optical  fibers  assume  the 
creation  of  an  electric-dipole  second-order  suscep¬ 
tibility  [1-5).  There  are  primarily  two  different 
mechanisms  proposed  to  explain  the  induced  'W’c- 
tric-dipole  susceptibility.  One  is  that  a  large  (  -  10~- 
105  v/m)  electric  field  is  induced  into  the  core  of 
the  fiber  producing  an  effective/'2’  from  an  allowed 
X(3)  interaction  (electric-field  induced  second-har¬ 
monic  generation).  The  second  mechanism  is  based 
on  aligning  defects  in  the  glass  through  the  intense 
polarized  laser  light  itself.  In  view  of  our  findings,  we 
argue  against  both  of  these  mechanisms. 

For  the  electric-field  induced  x<:>,  it  is  reasonable 
to  assume  that  the  dc-field  is  only  in  one  direction 
(,v  for  example)  which  is  parallel  to  the  polarization 
of  the  light.  It  that  is  so,  =/,[;,’■  =/[,%!  =0.  This 
can  be  deduced  from  the  relationship  /[,;’  = 
J'/.llI EjEkEfc  [15],  Instead  (see  table  4)  we  ob¬ 
serve  that /[,2J  ~  0.2/.[.;(  and  -  0. 1  In  the  dc- 
field  picture,  this  means  that  light  polarized  in  the  .v- 
direction  should  generate  a  dc-field  of  -  1CP  V/m  in 
the  .v-direction  and  ~  103  V/m  in  the  v-direction  si¬ 
multaneously.  We  find  this  mechanism  highly  im¬ 
probable. 

We  acknowledge  the  recent  measurements  of 
Mizrahi  et  al.  [16]  where  they  find 
0.00S/[;J  in  good  agreement  with  a  dc-field  picture. 
There  is.  however,  one  major  difference  between 
these  two  experiments  that  we  believe  can  explain 
the  different  results. 

In  the  experiment  of  Mizrahi  et  al.,  external  seed¬ 
ing  at  the  second-harmonic  wavelength  was  used  with 
a  peak  power  of  20  mW  polarized  in  the  .v-direction. 
In  our  experiment,  we  used  internal  seeding  which 
had  a  peak  power  between  -  10  "6- 10"  *  W,  de¬ 
pending  on  which  of  the  two  fibers  that  was  used, 
and.  furthermore,  this  internally  generated  green  light 
had  a  third  its  power  in  the  y-direciion  when  ex¬ 
cited  in  the  .v-direction  (tabic  3). 
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If  we  now  assume  that  SH  in  optical  fibers  occurs 
through  the  same  physical  mechanism  whether  in¬ 
ternal  or  external  seeding  is  used  in  the  preparation 
stage,  it  is  dear  that  the  electric  field-induced  model 
cannot  support  both  these  experimental  observa¬ 
tions.  In  the  case  of  external  seeding,  it  seems  that 
using  very  strongly  polarised  green  light  in  the  prep¬ 
aration  stage  suppresses  the  inherent  possibility  for 
growth  of  the  perpendicularly  polarized  green  lighu- 
If  the  mechanism  of  aligning  defects  were  correct,  we 
would  again  not  expect  such  a  large  growth  oi  SH 
light  perpendicular  to  the  preparation  polarization 
as  we  observe. 

Based  on  our  measurements,  we  suggest  that  the 
intrinsic  electric  quadrupole  second-order  suscepti¬ 
bility  is  the  origin  for  the  photoinduced  SH  light  and 
that  growth  occurs  through  improved  phase-match¬ 
ing.  We  justify  this  through  the  remarkably  small 
change  in  the  and  xl-U /Xxll  ratios  com¬ 

pared  with  the  overall  growth  of  SH  light  ( 7— S  or¬ 
ders  of  magnitude).  We  further  postulate  that  the 
improved  phasematching  is  due  to  periodic  altera¬ 
tions  in  the  electric  quadrupole  susceptibility  along 
the  core-cladding  interface  so  as  to  cancel  the  phase- 
mismatch.  In  principal,  a  periodic  refractive  index 
couid  account  for  the  improved  phasematching.  In 
this  case,  the  Tj.y'  would  be  multiplied  by  a  factor  i, 

( Sn  k  A )  where  J,  is  a  first-order  Bessel  function.  i\n 
is  the  amplitude  of  the  periodic  change  in  refractive 
index,  k  is  proportional  to  r[2>  and  A  is  the  peri¬ 
odicity.  Tspica!  numbers  for  an  optical  fiber  give  that 
7,  is  four  orders  of  magnitude  too  small  for  this  to 
be  feasible.  If.  instead,  the  nonlinear  susceptibility  is 
modulated,  r,:>  will  be  multiplied  by  a  factor  close 
to  one  in  magnitude.  There  have  been  reports  of  pe¬ 
riodic  structures  in  prepared  fibers  [17],  but  no  ex¬ 
periment  has  been  performed  showing  what  physical 
parameter  is  modulated. 

In  support  of  our  “phase-matching  model",  we 
would  like  to  draw  attention  to  the  recent  discovery 
that  the  third-harmonic  signal  in  an  optical  fiber  can 
alsz>  grow  with  time  [18].  To  first  order  this  growth 
was  not  accompanied  with  an  increase  of  the  third- 


order  nonlinear  susceptibility.  This  we  interpret  as 
another  indication  that  the  growth  of  harmonic  sig¬ 
nals  in  optical  fibers  are  primarily  from  improved 
phase-matching. 

In  conclusion,  we  have  measured  the  ratios  of  in¬ 
trinsic  and  induced  nonlinear  susceptibilities  in  two 
different  optical  fibers  as  well  as  the  polarization  de¬ 
pendence  for  the  SH  light.  In  addition,  we  have  pro¬ 
posed  a  new  approach  in  explaining  how  SH  light 
can  grow  in  an  optical  fiber. 
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